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Abstract 

We compute euclidean coordinate space correlation functions in a holographic model of QCD. 
We concentrate, in particular, on channels that are related to the U{\)a problem, the flavor-singlet 
axialvector, pseudoscalar meson, and pseudoscalar glueball (topological charge) correlator. We 
find that even a very simple holographic model defined on a slice of AdS^ provides a qualitatively 
correct description of QCD correlation functions. We study the role of anomaly terms, and show 
that both euclidean positivity and low energy theorems based on the axial anomaly relation are 
correctly implemented. We compare the results with expectations from an instanton model of the 
QCD vacuum. 
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I. INTRODUCTION 



Recently, much progress has been achieved in reahzing the old idea that large gauge 
theories are related to string theory The crucial development in this regard was 

the discovery of the correspondence between A/" = 4 super Yang-Mills theory and type 
IIB string theory on AdS^ x 5*5 [3]. Current effort is aimed at finding string duals of 
theories more closely related to QCD, in particular theories that exhibit asymptotic freedom, 
confinement and chiral symmetry breaking. Two avenues of research are being pursued. The 
top-down approach is based on the study of generalizations of the original AdS / CFT setu 
that incorporate fundamental fermions, chiral symmetry breaking and confinement ^,[5, 
I?]. The bottom-up approach is founded on models of QCD that are defined on warped higher 
dimensional spaces and incorporate the general principles, in particular holography, of the 
AdS/CFT correspondence j^, y, Q. The aim of the bottom-up approach is to find model 
independent features of the holographic description, to provide guidance for the top-down 
approach, and to develop a geometric language for thinking about gauge theories like QCD. 

Most of the published works on holographic models focus on the hadronic spectrum or on 
hadronic form factors. In this work we shall study euclidean correlation functions in a sim- 
ple holographic model of QCD. Euclidean correlators provide a bridge between perturbation 
theory and the operator product expansion, which control the short distance behavior, and 
the mass gap, which governs the long distance behavior llj. Phenomenological informa- 
tion on euclidean correlation functions comes from lattice data and from spectral functions 



12 



3. 



measured in processes like e"'"e~ annihilation into hadrons or hadronic r decays 

We will focus, in particular, on correlation functions related to the U{1)a problem. 
This includes the flavor-singlet axialvector, pseudoscalar meson, and pseudoscalar glue- 



Dall channel. 

Q 
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'he correlators in these channels have been studied in the instanton model 
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171], and there are some recent results from lattice simulations 
There are a number of qualitative questions related to the U{1)a sector in QCD that are 
still not fully understood. These include the origin of the topological susceptibility in pure 
gauge QCD, the microscopic mechanism that causes the susceptibility to vanish in QCD 
with massless fermions, and the role of instantons or other topological objects. The present 
work is motivated by the idea that holographic models of QCD can shed some light on these 
questions. 
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II. HOLOGRAPHIC MODEL 



We consider the holographic model introduced by Erlich et al. [10| and extended to the 



flavor singlet sector by Katz and Schwartz 
lagrangian 



2l|. The model is defined by the 5-dimensional 



S = Jd''x^i^-^TT[Fl + F^^)+TT{\DX\' + 3\X\ 



(1) 



^\DY\^+'-^[Y''f det{X) + h.c. 
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where X = XH"' is a scalar field, F is a complex fiavor singlet scalar, F^y = d^A^ — 
dpA^ — i[A^,Ay] (for L/R) is the field strength tensor corresponding to the gauge field 
= ^LR^"'- Here, are the generators of U{3)f. The covariant derivative is -D^X = 
dfj,X — iAif^X + iXAr^. The model is defined on an AdS^ metric 

ds^ = ^ (-dz^ + dx^dx^ , (2) 

with a "hard wall" cutoff Zm (0 < z < z^). The 5-dimensional masses of the fields are 
determined by the correspondence between fields on AdS^^ and operators on the boundary, 
m\ = (A — p)(A + p — 4). Here, ms is the 5-d mass, and A is the dimension of a p-form 
operator on the boundary. The scalar field X corresponds to the operator qiqj with A = 3 
and the complex scalar field Y corresponds to g^iGG + iGG) with A = 4. The gauge field 
has p = 1 and couples to the A = 3 operator qi'^^qj. The terms in the action that only 
involve Fl^r and X preserve the full U (3)_2< U (3) symmetry [l^ , and the last term generates 



the anomalous Ward identities of QCD 211]. 



The expectation values of X and Y are determined by classical solutions to the equations 
of motion in the presence of sources M and c, where M is the quark mass matrix, and c is 
related to the strong coupling constant (see equ. ([9]) below). We have 

{X,j) = a,,z'' + Mi,z, (r) = S 2;^ + c, (3) 

where and S correspond to the vacuum expectation values (jiiqj) and (G^) of operators 
conjugate to the sources M and c. In a more complete, top-down, model the expectation 
™lu« . and S a.e deternJned dynanncaUy fl Q. Alternatively, one n.ay inco.po.ate the 
dynamics of chiral symmetry breaking into a boundary action for X and Y [9|. Here, rather 
than specify these boundary terms, we shall take a and H to be free parameters of the 




FIG. 1: Vector and axialvector current correlation functions. We show the ratio of the correlation 
function to the free correlator as a function of euclidean separation r. The solid curves show the 
result in the holographic model. The data points are taken from an analysis of Aleph data on 



hadronic tau decays 



1^. 



model. In this work we are not specifically interested in flavor mixing (see 2l|) and we will 



set cTjj = 5ija and Mij = dijUi. 



The coupling constants and c are determined by matching the short distance behavior 
of correlation functions to QCD. Consider the vector current correlation function 



where = q'j^t°'q and = —q^. Using the AdS/CFT dictionary we get 



V{q,z)d,V{q,z) 



(4) 



(5) 



where V{q, z) is the bulk-to-boundary propagator, V^{q, z) = V{q, z)Vl^{z), and V^(g) is the 
boundary value of the field. The propagator has to satisfy the equation of motion. In the 
Vz{x, z) = gauge the linearized equation of motion for the gauge field is 



dJ-dMq,z))+-V^{q,z) 



0. 



(6) 



,;z J z 

A solution on the complete AdS^ space with z — > 0) = 1 is given by z^ = 
— ^{qz)Yi{qz) and 

Ily{Q') = -^\og(Q'). (7) 



This result can be compared to the perturbative one-loop contribution ny(Q^) = 
—Nc/{24:T{'^) ■ log((5^). This leads to the matching condition 

Ol = (8) 

We observe that ~ 1/iVc which shows that vector mesons at large Nc are weakly coupled, 
as expected from general large Nc arguments. The constant c is fixed by matching the fiavor- 
singlet axialvector correlation function and using the axial anomaly relation. The result is 



m 



a. 



and we shall follow 



2A'/^. (9) 



21| and allow to run as a function of z according to the one-loop 
beta function. Having fixed the coupling constants we can now compute the correlation 
functions. For this purpose we shall employ the representation of the Green function in terms 
of eigenfunctions of the five-dimensional Sturm-Liouville problem. In the vector channel 



nv(Q^) = E . /i'' . X (10) 



T ml, (g2 + 



where = m?p f^ is an eigenvalue of equ. ([6]) subject to the boundary conditions V^(0) = 
and V'{z„i) = 0. The eigenfunctions are normalized according to 

dz^Vk{z)Vi{z)=6,i, (11) 



and the decay constants fp^k are given by 



(12) 



The eigenmode representation ffTOl) also determines the spectral function and the euclidean 
correlator. The spectral function is a sum of delta-functions, and the euclidean correlation 
function is 

ny(r) = Y.fl,D{mp,k,r), D{m,T) = K,{mT), (13) 

k 

where D{m,T) is the euclidean propagator of a scalar particle with mass m. In the vector 
channel the spectrum can be obtained almost in closed form. The eigenvalues rup^k are 
determined by the zeros of the Bessel function of the first kind and order zero, Jo("^p,A:^m) = 
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FIG. 2: Pseudoscalar octet and singlet current correlation functions. The correlation functions are 
normalized to the correlator in the free theory. 

0. Using the location of the first zero and the physical mass of the rho meson we can set 
the scale, = 0.62 fm. The coupling constant of the k'th resonance is given by 

/p,fc = Y( V- (1'^) 

The euclidean correlation function is shown in Fig. [H We plot the ratio of the euclidean 
correlator T{y{T) over the correlator in the non-interacting theory, liyij) = Nc/{3tt'^t^). 
Because of asymptotic freedom this ratio approaches one as r —>■ 0. The result can be 
compared to the experimental results compiled in 1^. The agreement is quite good, but 
at short distances perturbative logarithms and power corrections are missing. The short 
distance behavior can be improved by adding extra terms to the action given in equ. ([1]), 
but this is beyond the scope of the present work. The agreement with experimental data is 
comparable to resonance saturation models constrained by Weinberg sum rules and other 



QCD inputs [221, l23|]. This is not surprising, as equ. (fTOll is equivalent to the resonance 
saturation ansatz, equ. (IHl) is the statement of local duality, and Weinberg sum rules are 
realized as a consequence of chiral symmetry. 
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III. PSEUDOSCALAR AND AXIALVECTOR CORRELATION FUNCTIONS 



The equations of motion for the transverse components of the axialvector gauge field is 

where we have introduced the notation v = mz + az^ for the vacuum expectation value 
of the scalar field. This equation differs from the equation for the vector field by an extra 
contribution due to (explicit and spontaneous) chiral symmetry breaking. The longitudinal 
component of the axialvector current mixes with the pseudoscalar current, and pseudoscalar 
fluctuations have be to taken into account. We define octet r/" (a = 1, . . . , 8) and singlet 
?7°, a fields by 

= {Xij) exp(^r/'^t'^), Y = (F) exp(2a). (16) 

We also define the longitudinal component of the axialvector field d^i^"" = — A^j^ (in 
= gauge). We first consider the flavor non-singlet sector. The equations of motion for 
the longitudinal axialvector yj" and the pseudoscalar field vr" are 

9,(^9,^'^) +^(77'^-^'^) = 0. (17) 



q'^d,^" + ^d^Ti'^ = 0. (18) 



z^ 



The solutions are normalized according to 



dz 4 < « - cpt) = 5,1 (19) 
z-^ 



and the pion decay constant is given by 



J-K,k 2 



(20) 

We consider Nc = ?) and set m = 2.23 MeV and a = (323 MeV)^ in order to reproduce the 
experimental values of and for the lowest pion excitation. The trace of the axialvector 
correlation function, 11^ = ^Afi^ is shown in Fig. [TJ The splitting between the vector and 
axialvector correlators is related to the difference in mass and coupling between the vector 
and axialvector mesons, and to the pion contribution. The latter dominates at large distance 
and causes 11^ to become negative. In the model considered here chiral symmetry is restored 
high in the spectrum (m^ m^^ as A; — oo) and the vector and axialvector correlators 
are very nearly degenerate for r < 0.5 fm. 



In the singlet sector there is also mixing with the pseudoscalar glueball field a. The 
equations of motion are 

,2 „2 



/„2 \ „2 Nf 

qh^.p" - gyaM' - 9lc'd.a = , 
where we have defined k = c^f kq. The normalization condition is 



dz 



Ski ■ 



(21) 
(22) 
(23) 

(24) 



For a normalized eigenmode the coupling constants to the axialvector, pseudoscalar quark- 
anti-quark, and pseudoscalar glueball current are given by 



Jri',k 2 



\',k — ^ 



dzVk 



hri',k 



95 z 

In Fig. [2] we show the pseudoscalar correlation functions 

Up{Q')6'''= I d'xe''^qe^,q{x)qt'^5qm 



2Nf 



(25) 



(26) 



in the fiavor singlet and non-singlet sector. In the non-singlet sector the correlation function 
is dominated by the light pion excitation. In the singlet sector the would-be Goldstone boson 
receives an extra contribution to its mass from the anomaly term. We have set k, = 20 and 
find rrir]' = 650 MeV. (For k, ^ 1 the mass of the rj' becomes weakly dependent on the 
value of K. We also note that the experimental result for the rj' mass, m^/ = 957 MeV, can 
be reproduced if the strange quark mass is taken into account.) Like chiral symmetry, the 
axial U{1)a symmetry is effectively restored in the highly excited part of the spectrum. For 
k ^ 1 the equations of motion for rj^, ip^ and a effectively decouple, and the solutions to 
Eqs. (12111231) alternate between solutions to the equations of motion in the non-singlet and 
the glueball sector of the theory. 

IV. TOPOLOGICAL CHARGE CORRELATOR 



The pseudoscalar glueball correlation function is 

np(Q') = - 



^4^gi<?x /^^c'(a;)^GG(0) 

OTT OTT 



(27) 
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FIG. 3: Pseudoscalar glueball correlation functions. The correlation functions are normalized to 
free field behavior. The curve labeled k = corresponds to the pure gauge theory without anomaly 
term, whereas the curve labeled k ^ includes the effects of an anomaly term. The curve labeled 



Note that as/ (Svr) ■ GG is the topological charge density. In order to study the topological 
charge correlator in the holographic model we begin with the equations of motion in the 
pure gauge sector of the theory. We have 



where k = nv'^f . The solutions of this equation are particularly simple in the limit c = const 
and K ^ 0. In this case the eigenvalues in the pseudoscalar glueball channel are given by 
the solutions of JiijnG^kZm) = and the coupling constants are 



The groundstate in the pseudoscalar glueball channel is heavier than the ground state rho 
meson by a factor ~ 1.5. The glueball correlation function is shown in Fig. [31 The shape of 
the correlation function is similar to the vector meson correlator, but since the mass scale 
is larger the correlator is smaller. The topological susceptibility is defined as the integral of 
the pseudoscalar glueball correlation function 



Nf = 3 shows the result in full QCD, including mixing with the rj' . 




(28) 




(29) 




(30) 
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FIG. 4: Integrand for topological susceptibility in pure gauge QCD without an anomaly term. The 
integrand changes sign at z ~ 1 fm, and the topological susceptibility vanishes. 

where Qtop is the topological charge and V is the volume. This integral diverges at short 
distance, and care has to be taken in defining the subtraction scheme. In the present case 
it is sufficient to subtract the contribution from the free gluon bubble, Hp = 6Q;s/(vr^x^). 
This corresponds to subtracting the correlator in the infinite (non-cutoff) AdS^ space. The 
integrand 27i'^x'^{Ilp{x) —11^(0;)) is shown in Fig. HI We observe that the integrand changes 
sign, and we find that the topological susceptibility vanishes. This result, as well as the 
subtraction scheme, can be checked by computing the topological susceptibility directly from 
the solution of the equation of motion (l28l) for = 0. The only solution that satisfies the 
boundary conditions a(0) = 1 and dza{zm) = is a{z) = 1, and the topological susceptibility 



^2 







(31) 

vanishes. This result does not depend on the simplifying assumption c = const. It does, 
however, depend on the choice k = 0. The correlation function for k 7^ is also shown 
in Fig. [3l The anomalous term suppresses the correlator and leads to a non-vanishing 
(and positive) topological susceptibility. In the simplest model, defined by equ. fl25]) . the 
topological susceptibility has a non-perturbative ultra-violet divergence. This can be seen 
by treating the anomaly term as a small correction. In this limit the topological charge 



correlator is 



2l| 



n-«' = -2iv7/o 



2 



(32) 
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FIG. 5: Five dimensional = modes in the pseudoscalar sector. We show solutions of the 
equations of motion for a{z) (blue, solutions start at a(0) = 1) and ri^{z) (red, satisfy r/'^(0) = 0) 
for different values of the quark mass m. The topological susceptibility is related to the negative 
of the curvature of a{z) near z = 0. 

where the expression in the square brackets is the bulk-to-boundary propagator of the a 
field in the infinite AdS geometry (for R = 0). For Q'^ = the integrand is singular near 
z = unless R,{z) vanishes faster than z"^ as z ^ 0. This problem can be understood by 
comparing equ. fl5^ with the expectation from an instanton model, see appendix \M In 
terms of instantons the integral dz / z^ arises from the integral over the size of the instanton, 
and the power z^ is fixed by classical scale invariance. The Bessel function corresponds to 
the Fourier transform of the topological charge density, and the order of the Bessel function 
is determined by classical scale invariance and the dimension of the operator GG. The 
instanton contribution in QCD contains an extra factor exp{—87i'^/g'^{z)), where Sn'^/g'^ is 
the instanton action, and g{z) is the running coupling constant. Asymptotic freedom implies 
that this term scales as (for z 0), where b is the first coefficient of the beta function. 
As a consequence, the integral is well behaved for small z, but diverges at large scale sizes. 
In the holographic model the infrared problem is cured by the hard wall cutoff, but the 
simplest version of the model has an ultraviolet divergence instead. 

The problem can be addressed by considering a more complicated functional form for 
R,{z). Here, we proceed directly to the unquenched theory, which has a finite topological 
susceptibility. The topological charge correlator in the theory with three fiavors is shown in 
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Fig. ([3]). There is an extra contribution which is dominated by mixing with the lowest rj 
resonance. This term increases the topological charge correlator at long distance, and tends 
to cancel the topological susceptibility generated in the pure gauge theory. The topological 
susceptibility can be determined most accurately by solving the equations of motion for 

(^^^«) + - = 0' (33) 
v^aV + c^^^a = 0, (34) 

subject to the boundary conditions a(0) = 1, ?7°(0) = 0, dza{z,m) = dzrp{zm) = 0. Solutions 
to equs. (13311341) for different values of the quark masses are shown in Fig. [51 We observe 
that the topological susceptibility scales with the quark mass, as expected from low energy 



theorems based on the axial anomaly 24, 25|. For m — 2.2 MeV we obtain Xtop — (70 MeV)3, 



in agreement with the expected value Xtop = ma/Nj. 



V. SUMMARY 



We have studied euclidean coordinate space correlation functions in a holographic model 
of QCD. The correlation functions in the flavor non-singlet vector and axialvector channel 
agree well (within ~ 20%) with experimental data from hadronic tau decays. The holo- 
graphic model does not describe perturbative logarithms and power corrections. In principle 
this can be addressed by modifying the AdS^ geometry, or by adding higher dimensional 
terms in the five-dimensional action. 

In the pure gauge sector, without an anomaly term, the model generates a non-trivial 
pseudoscalar glueball spectrum but the topological susceptibility is zero. The anomaly term 
gives a negative (repulsive) contribution to the pseudoscalar glueball correlator and leads 
to a non-zero topological susceptibility. In the full theory an extra attractive contribution 
arises from mixing between the pseudoscalar glueball and the pseudoscalar meson (?]') field. 
These two terms cancel as the quark mass goes to zero, and the topological susceptibility 
vanishes. The attractive term is longer ranged than the repulsive one, so the vanishing 
of Xtop can be viewed as being due to topological charge screening. It also interesting to 
compare the results to expectations from the instanton model. The anomaly contribution to 
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both the topological charge and the pseudoscalar meson correlator agree with the structure 
of the dilute instanton gas result. 

In this work we studied a very schematic model based on a slice of AdS^. Clearly, it 
is of interest to study correlation functions in "top-down" models. A very detailed study 



of anomalous correlators in the A/" = 4 theory was carried out by Dorey et al. 26| Some 
attempts to study topology and the y' meson in extensions of the simplest AdS/CFT setup 



32l |. It is also important to further clarify the role of 



can be found in 



21, 



28|, 



29 



31 



instantons in the large limit. Witten argued that instanton effects are suppressed in 
the weak coupling limit of a large field theory 33|]. In QCD the situation is not clear, 
because the theory is classically conformal (it has instantons of all sizes), and in the large Nc 
limit only small instantons are suppressed. We have previously argued that the instanton 



size distribution in the large Nc limit might be a delta function [3J]. In the language of 
holographic QCD this corresponds to an anomaly term k,{z) which is sharply peaked at 
some critical distance z* in the fifth dimension. 
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APPENDIX A: CORRELATION FUNCTIONS IN THE FIELD OF AN INSTAN- 
TON 



In this appendix we collect a few results for the single instanton contribution to hadronic 
correlation functions. The instanton contribution to the pseudoscalar glueball correlation 



function is 



15 



16| 



p5 



n 2 



Q'p'K2{pQ) 



(Al) 



np(Q) = - 

where the Bessel function K2{pQ) arises from the Fourier transform of the topological charge 
density of an instanton. The factor 2 arises from adding the contribution of instantons and 
anti- instantons. The instanton size distribution is given by 35 1 



d{p) 



C 



2Nc 



Nc 



exp 



n ^}p^ 

f 



(A2) 



13 



47r2 

m*fp = mp- —{qq)p\ 



C, 



Nc 



3 

0.466 exp(-1.679A/'e)1.34^/ 



(iV,-l)!(iV,-2)! ' 
where g{p) is the one- loop running couphng constant 

4fr = -Mog(pA), b=^N,- ^Nf. (A3) 
9 [P) 3 3 

The factor 1/p^ is related to classical scale invariance and matches the determinant of the 
metric tensor on AdS^. We also note that the factor {m* p)^^ matches the factor v^^ = 
{mz + az^)'^f in equ. (155]) . The remaining terms in d{p) are related to fluctuations around 
the classical instanton solution and have no obvious counterpart in the holographic model. 

The instanton contribution to the flavor non-singlet pseudoscalar meson correlation func- 
tion is [36 1 

This result includes only the contribution of the fermion zero mode in the quark propagator. 
The instanton contribution to the flavor singlet correlator is the same up to an overall sign. 
The Bessel function Ki{pQ) is the Fourier transform of the fermion zero mode. The order is 
related to the dimension of the Fermion field, and matches the bulk-to-boundary propagator 
for the pseudoscalar field on AdS^. 

The instanton contribution to the vector current correlation function is 37 1 

n,.(Q) = {q,q. - 6,.q') ^ J ^ d{p) |^ " ^ J^x K, (2pQ/(l - x')'/')^ 



n.(Q) = 2 / ^ dip) -1 [QpK,ipQ)f . (A4) 



(A5) 



The first term is a pure power correction which corresponds to the {g^G"^) and m{qq) terms 
in the OPE. These terms are not present in the holographic model. The second term, 
which is due to non-zero modes of the fermion propagator in the instanton field, also has no 
counterpart in the model considered here. 
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